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ABSTRACT 


An analytical evaluation was made of five turbine blade cooling configurations for 
advanced air breathing engines on the basis of the time to a local creep strain of 1 per- 
cent for cast IN 100. The longest life was shown by the transpiration design (neglecting 
oxidation problems) followed by the multiple small hole design. Curves were established 
to predict the multiple small hole blade life under various combinations of turbine inlet 
and coolant temperatures and coolant to total gas flow ratios. 



STEADY -STATE STRESS RELAXATION ANALYSIS OF TURBINE 
BLADE COOLING DESIGNS 
by Albert Kaufman 
Lewis Research Center 

SUMMARY 

An analytical evaluation was made of five turbine blade cooling configurations for 
advanced airbreathing engines on the basis of the steady-state operating times to 1 per- 
cent creep strain. The configurations consisted of a strut insert blade with horizontal 
fins, a film and convection cooled, a transpiration cooled, a multiple small hole, and a 
liquid metal cooled design. Analytical temperature distributions for a hypothetical 
Mach 3 supersonic cruise aircraft engine were used. The blades were assumed to oper- 
ate under the centrifugal loading at cruise conditions of a typical first stage turbine ro- 
tor blade. Blade lives and stress and strain distributions were calculated by means of a 
stress relaxation computer program which used creep data for cast IN 100, The tran- 
spiration cooled design showed the longest operating life followed by the multiple small 
hole design. However, the latter was concluded to be more practicable for operation at 
high metal temperatures because of oxidation problems with the porous shell used in 
transpiration cooling. The superiority of the multiple small hole over the strut insert 
and film and convection cooled blades in the analysis was believed due largely to a more 
nearly optimum distribution of the coolant flow over the blade. 

Design curves were established for the multiple small hole design to predict lives 
to 1 percent creep strain for various combinations of turbine inlet and coolant tempera- 
tures and coolant flow to total air flow ratios. 

INTRODUCTION 

Advanced air breathing engines with high turbine gas inlet temperatures require 
cooling of the turbine blades. In this study, a number of cooling schemes have been 


evalxxated on the basis of the steady-state creep lives under supersonic cruise operating 
conditions. 

Present blading materials are limited to metal temperatures of approximately 
1800° F (1255 K) because of strength and oxidation considerations. However, for some 
engines the compressor discharge air which is the most convenient source for a cool- 
ant is already at a temperature above 1000° F (811 K). In order to utilize this coolant 
with a minimum amount of compressor bleedoff and cool the blade to within the limits 
on temperature levels and gradients set by the applied loads and the material strength 
and oxidation properties requires cooling configurations of high effectiveness. 

A number of cooling schemes were analyzed from a heat transfer point of view in 
references 1 to 3. The engine operating conditions for this analysis were Mach 3, 

75 000 feet (22. 86 km), and a 2270° F (1517 K) gas turbine inlet temperature with a 
1270° F (961 K) blade coolant inlet temperature. The blades had a chord width of about 
2 inches (5. 1 cm), a span length of 5. 3 inches (13. 46 cm), a hub to tip diameter ratio of 
0. 81, and were assumed to be fabricated from cast IN 100. Temperature distributions 
were determined analytically for a number of airfoil sections at various locations along 
the blade span; these results are reported in references 1 to 3. These blades were de- 
signed to sustain 1000 hours of steady-state operation based on three standard deviations 
of stress rupture properties and an empirical design procedure. In reference 4 cyclic 
lives based on a proprietary analysis were presented for sea level transient conditions. 

The purposes of this study are threefold: (1) to present a method for examining 
various types of cooled blades under steady-state conditions taking into account the 
stress relaxation due to creep; (2) to investigate the steady-state lives of a number of 
different cooling configurations from references 1 to 3 using a stress relaxation method; 
and (3) to establish preliminary design curves for potential lives of the more promising 
of these configurations under various combinations of coolant flow ratio, coolant temper- 
ature, turbine inlet gas tender ature, and centrifugal force. 

In order to evaluate these designs, the centrifugal loads were based on operating 
conditions for a first stage rotor blade having a tip speed of 1243 feet per second (383 
m/sec) and a hub-tip diameter ratio of 0. 75. The steady-state operatii^ lives to 1 per- 
cent creep strain were determined for five of the blade cooling conf^rations analyzed 
in references 1 to 3. These five schemes were (1) a strut inserted design where the mid- 
chord region is convectively cooled through horizontal fins, the leading edge is impinge- 
ment cooled, and the coolant is discharged through a split trailing edge, (2) a film 
convection-cooled design where the midchord region is convection cooled and the leading 
and trailing edges are convection and film cooled, (3) a transpiration cooled design using 
a strut supported porous shell, (4) a multiple small hole design where the cooling is done 
primarily by film cooling by ejection of the coolant from the airfoil surface, and (5) a 
liquid metal cooled design. Temperature distributions from references 1 to 3 for an ap- 
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proximately 0. 026 coolant to total gas flow ratio were used for this analysis. 

The creep lives of these blade cooling conf^rations were calculated by means of a 
computer program which has the capability of taking into account the mechanical 
stresses from the external centrifugal and gas loads, the thermal stresses due to the 
temperature gradients, and the redistribution of stresses with time because of creep. 

One effect of creep is to cause a relaxation of the high thermal stresses. If this relaxa- 
tion is ignored the predicted blade life will be too pessimistic. On the other hand, if the 
thermal stresses are ignored entirely the predicted blade life will be too optimistic. 
Therefore, a stress relaxation analysis is required to realistically evaluate blade life. 

The analysis does not consider effects of damage due to thermal fatigue, vibrations, 
impact, or erosion. Some or all of these effects can occur at various times durir^ the 
service lives of turbine blades under engine design, takeoff, and off-design conditions. 
The actual blade lives will depend on the Interaction of all these factors with the steady- 
state conditions. Even if the magnitudes of all these nonsteady- state effects were known, 
a prediction of blade lives combinii^ these effects with the steady-state effects is be- 
yond the current state of the art. 

Under steady-state conditions, the elapsed time to obtain 1 percent creep strain 
(henceforth called creep life), the Initial and final stress distributions, and the final 
strain distributions are presented for each design. Creep lives for these designs are 
compared on the basis of the initially (prior to stress relaxation) worst combination of 
stress and temperature, the centrifugal stress and the average temperature, and the ef- 
fect of creep as determined from the stress relaxation computer program. 

Conclusions are drawn as to the high temperature capabilities of these designs. De- 
sign curves are established in order to evaluate the effect on the creep life of the mul- 
tiple small hole configuration of various combinations of coolant temperature, turbine 
inlet gas temperature, cooling air flow to total gas flow ratio, and centrifugal force. 


SYMBOLS 


A 

E 


^3 


F 

k 

L 

M 


. 2 2 
area, m, ; m 

coefficients defining a plane 

2 

modulus of elasticity, psi; N/m 
centrifugal force, lb; N 
curve fix constants 
span lei^th, in. ; m 

bending moment with respect to original axis, lb- in. ; N-m 
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m’ bending moment with respect to centroidal axis, lb- in. ; N-m 

n number of nodes 

R radius of curvature, in, ; m 

T temperature above 70° F, °F; K 

Tjj temperature, °R; K 

t time, hr 

w deflection, in. ; m 

X distance along x axis, in. ; m 

X distance to centroid in direction of x axis, in. ; m 

X’ distance aloi^ x' centroidal axis, in. ; m 

y distance along y axis, in. ; m 

y distance to centroid in direction of y axis, in. ; m 

y' distance along y* centroidal axis, in. ; m 

z distance along z axis, in. ; m 

a coefficient of thermal expansion, in. /in. /°F; m/m/K 

e strain, in. /in. ; m/m 

0 angle between x’ and x axes, rad 

2 

a stress, psi; N/m 

Subscripts: 
c creep 

k node index 

0 hub 

p plastic 

t total 

X, y, z with respect to x, y, z axes, respectively 
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ANALYTICAL CONDITIONS 


Cooling Configurations 

The five cooling schemes that were analyzed are illustrated in figures 1 to 5. hi the 
strut insert des^n of figure 1, air flows up the central cavity formed by the strut insert 
and through holes at the leadii^ edge of the insert to impingement cool the blade leading 
edge. The air then circulates through horizontal fins between the shell and strut and 
discharges through slots in the trailing edge. For the purpose of this analysis, the 
outer shell, strut, and horizontal fins are assumed to be integral. 

The film and convection cooled design of figure 2 has the cooling air injected through 
the blade base into two central and one leading edge cavities. The air then circulates up 
and down a series of vertical passages. At the leading edge the air passes through 
small holes in the wall of the adjacent vertical passage, impinges on the internal leading 
edge surface, and discharges through film cooling holes. The trailing edge is convec- 
tively cooled by air discharging throi^h slots. 

The transpiration cooled des^n of figure 3 has a shell of wire-form porous material 
attached to a strut which carries the shell. Air flows up the central plenum of the hol- 
low strut, through metering holes to the strut surface which consists of a waffle-like 
pattern of cavities, and then through the porous shell. The shell material is cooled by 
a combination of convection and film cooling which is highly effective because of the 
many fine pores through the shell structure. 

The multiple small hole des^n of figure 4 is essentially a film cooling scheme which 
employs a large number of closely spaced holes in the shell wall. These holes are con- 
siderably larger than the pores in the fine wire mesh used in the transpiration cooled de- 
sign of figure 3 and are, therefore, less susceptible to clogging because of metal oxida- 
tion at high temperatures, hi the multiple small hole des^n the shell is stiffened by 
cross ribs and is capable of supporting itself imder engine operatiig conditions without 
a strut. 

In the design shown in figure 5, the leadiig e<%e and midchord regions are convec- 
tively cooled by employing a liquid metal; in this case, potassium. Boiling of the po- 
tassium occurs as it flows through radial passages inside the airfoil part of the blade. 
Condensation occurs below the blade platform in heat exchangers which are cooled by 
compressor discharge air. The trailing edge is convectively air cooled. 
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^-68 slots, 0.03 by 0.06 in. (0. 76 
by 1.52 mm), 0.02 in. (0.51 mm) 
between slots, suction side insert 






Figure 2. - Film -convection cooled design. 
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-ir 


T 


■0.012 in. 
(0.305 mm) 
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p 
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!° 

/ o 

4.004 in. 

{0.102 mm) diam 


uonriguration K 


/o 


Q "^^"0.008 in. (0.203 mm) 
p ^ ~~{~0 .016 in. (0.406 mm) 


/ o 

H ho. 


008 in. (0.203 mm) 


-0.004 in. (0.102 mm) diam 


Configuration C 


O- 

O 

O 

O 


L 

0.012 in. (0.305 mm) 
To. 024 in. (0.61 mm) 



0.008 in. (0.203 mm) 


^0.008 in. (0.203 mm) diam 
tapers from 0.004 in. 

(0. 182 mm) at insidewall 


0.010 in. (0.254 mm) diam holes 

spaced at 0.020-in. (0.51mm) CD-10206-33 


Figure 4. - IViultiple small-hole design. 
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Temperatu res 


The heat transfer analysis of references 1 to 3 was done for operating conditions as 
found in an advanced airbreathing engine with a 2270° F (1517 K) turbine inlet tempera- 
ture at Mach 3 and 75 000 feet (22 859 km) altitude. The blade material was assumed 
cast IN 100, The heat transfer calculations of references 1 to 3 cover a wide range of 
coolant flow ratios. In order to arrive at a uniform basis for comparing the designs, 
coolant flow ratios of approximately 0, 026 were selected, hi the case of the multiple 
small hole design where no results were presented for a 0, 026 coolant flow ratio, the 
ten^erature distributions were obtained by averaging the results for 0. 0.286 and 0, 0223 
coolant flow ratios, 

A survey was made to determine the critical span location for each design by calcu- 
lating a nominal stress rupture life based on the average temperature at the section and 
the centrifugal stress distribution. The average temperatures were obtained from the 
temperature distributions which were given for each configuration at various span loca- 
tions in references 1 to 3, It was found that in all cases the hub was critical when the 
temperature distribution at the section closest to the hub section of the airfoil was con- 
sidered to apply at the hub also. These sections were 8 to 20 percent of the span dis- 
tance from the airfoil base. The only exception to this was in the case of the transpira- 
tion cooled design where strut temperatures were only presented for the 60 percent span 
section. Since the strut is the load carrying member and the shell temperatures showed 
little variation with span distance, the 60 percent span temperature distribution was also 
applied to the hub section of the strut. This was a conservative approach since the hub 
section of the strut would be cooler than the 60 percent span section. The complete tem- 
perature distributions were given in references 1 to 3 for one section only, usually in the 
vicinity of the midspan. At other sections, metal temperatures were given at approxi- 
mately 15 points. Temperatures at intermediate locations for the purpose of this anal- 
ysis were obtained by a combination of interpolation and the use of the midspan temper- 
ature distribution as a guide. The ten^eratures that were used from references 1 to 3 
are shown in figures 6 to 10 for each design. 
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CD-10375-01 

Figure 9. - Temperature distribution for multiple small-hole design. 
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Figure 10. - Temperature distribution for liquid-metal- and convection -cooled design. 


Applied Loads 

For purposes of evaluating these configurations the centrifugal loads were based on 
the operating conditions and the geometry of a first stage turbine blade with a tip speed 
of 1243 feet per second (383 m/sec) and a hub-tip diameter ratio of 0.75. From analy- 
ses of similar configurations for the same turbine rotor blades, average centrifugal 
stresses of 31 000, 29 000, 37 000, and 30 000 psi (21 374, 19 995, 25 511, and 20 684 

9 

N/cm ) were assigned to the strut insert, film and convection cooled, transpiration 
cooled, multiple small hole, and liquid metal and convection cooled designs, respec- 
tively. These stresses were applied to the airfoil sections closest to the airfoil base, 
hereafter called hub sections. The effect of bending moments due to the gas forces was 
neglected in this analysis because it was assumed these would be cancelled out by tilting 
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the blades. However, bending moments due to the offsetting of the line of action of the 
centrifugal force as a result of blade bendir^ due to creep were considered. 


ANALYTICAL PROCEDURE 

Stress distributions in the airfoil sections were calculated by a stress relaxation 
procedure which is explained in detail in appendix A. This general program takes into 
account the applied loadii^ on the blade, thermal stresses, and creep effects under 
steady-state operation, hi accordance with present practice beam theory is used; that 
is, plane sections are assumed to remain plane. The program has the capability of 
treatii^ Instantaneous plastic flow or residual stress effects. However, all the stresses 
were in the elastic rar^e for all the cases considered in this investigation. 

In appendix B a procedure is suggested for determining the effect of creep on the 
bending moments. The effect of these moments was approximated by (1) assuming the 
deflection distribution over the blade span as a function of the slope (Cg or Cg in 
eq. (Al)) of the hub section on the assumption of small deformation, (2) determining the 
centrifugal force distribution over the span for a normal tapered airfoil shell, and (3) 
calculating the total restoring moment on the hub section as a function of the slope by 
multiplying (1) and (2) and integratiig the results. 

The analysis does not consider the effects of nonlinear spanwise temperature grad- 
ients. hi the cases of the transpiration cooled blade where the strut was assumed to be 
at a constant temperature and the multiple small hole blade where the spanwise variation 
of average section temperatures was fairly linear the analysis of appendix A is reason- 
ably accurate. The other configurations also had fairly linear spanwise temperature 
gradients in references 1 to 3 except near the tip where the temperatures leveled off or 
decreased rather sharply. The effects these have on the accuracies of the calculations 
are unknown. 

The basic nodal breakup of the airfoil sections for each configuration as shown in 
references 1 to 3 was used for this investigation. However, the strut for the transpira- 
tion cooled design had to be divided into small area elements since no nodal breakup was 
given for it in references 1 to 3. Distances from the centroids of each of these nodes to 
an arbitrary set of coordinate axes were determined. Temperatures based on the heat 
transfer analyses were assigned to each node. The material properties are tabulated in 
the program as functions of temperature and are, therefore, known at each node. 

Creep data were obtained for IN 100 alloy from the International Nickel Company for 
creep strains of 0. 1, 0. 2, 0. 5, and 1.0 percent at ten^eratures of 1350°, 1500°, 1700°, 
1800°, and 1900° F (1005, 1089, 1200, 1255, and 1311 K) and times from 10 to 1000 
hours. These creep data were correlated using a Larson- Miller parameter as shown in 
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Larson-Miller parameter, Tp (25.56 + logjQt) 
Figure 11. - Creep data correlation for IN 100. 


figure 11. The parametric constant was optimized from the available data by using the 
computer program described in reference 3. It was assumed that compressive stresses 
produce compressive creep strains of the same magnitude as the tensile stresses for 
any given absolute stress, time, and temperature. An equation of the form 

^°Sl0^c " ^1^ ^ [^2^ ^3] ^ *^4 

was fitted to the creep curves of figure 11. 

A strain hardening rule was used to determine cumulative creep strains. This 
procedure and the application of the rule are explained in appendix A. 



RESULTS AND DISCUSSION 


Creep Lives 

The creep lives to 1 percent creep strain are shown in table I for the Mach 3, 

75 000 feet (22. 859 km), 2270° F (1517 K) turbine inlet temperature, and 1270° F 
(961 K) blade coolant inlet temperature conditions. In addition to the results from the 
stress relaxation analysis shown in the third and fourth columns, creep lives are given 
for the worst combination of stress and temperature using initial conditions (0 hr at the 
steady-state conditions) in the second column and for the average temperature and stress 
at the hub section for each blade design in the last column. Creep lives from the stress- 
relaxation analysis are presented with and without the effect of restoring moments due 
to centrifugal force beii^ considered. 

Since the creep correlation curves of figure 11 were based on -data from 10 to 
1000 hours, results outside this range and, particularly, beyond 10 000 hours can be 
considered as only qualitatively significant. More accurate creep lives could have been 
obtained if the Larson- Miller correlation had covered creep data in the 1000- to 
100 000-hour range. Because of this deficiency in the creep data correlation, the 
results shown in table I are used only to compare designs and describe trends. 

The blade creep lives of various designs based on the stress relaxation analysis 
with restoring moments being considered were, in descending order: (1) the transpira- 
tion cooled design with an infinite life; (2) the multiple s mall hole design with 1 1 600 
hours; (3) the strut insert des^n with 7200 hours; (4) the film and convection cooled de- 

TABLE I. - RESULTS OF CREEP LIFE ANALYSES 


[Cruise speed, Mach 3; cruise altitude, 75 000 ft (22, 859 km); turbine inlet temperature, 
2270° F (1517 K); blade inlet coolant temperature, 1270° F (961 K), ] 


Blade cooling design 


Time to 1 percent creep strain, hr 



Based on initial 

Based on stress relaxation analysis 

Based on average 


conditions 

Uncorrected for 

Corrected for 

conditions 



restoring moment 

restoring moment 




due to centrifugal 

due to centrifugal 




force 

force 


Strut insert 

2430 

4 700 

7 200 

47 900 

Film convection 

186 

1 750 

2 850 

46 700 

Transpiration 

2530 

47 900 

Infinite 

Infinite 

Multiple small hole 

4800 

8 600 

11 600 

33 500 

Liquid- metal convection 

68 

130 

134 

498 

1 
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sign with 2850 hours; and (5) the liquid metal and convection cooled design with 134 
hours , 

The effect of the restoring moments was to cause an increase in the creep lives. 

The reason for this is that it counters the tendency of the blade section to rotate due to 
creep. The life of the transpiration cooled design increased from 47 900 hours to an in- 
finite life while that of the liquid metal design increased only slightly due to the restor- 
ing moments. However, the qualitative ranking of the blade designs in terms of creep 
life did not change. 

The reason for the superiority of the transpiration cooled design was that the strut 
had a relatively low average temperature of about 1400° F (1033 K) which more than 
compensated for the high centrifugal stress level required to support the porous shell. 
Similarly, the liquid metal design had such low life because the hub section was at an 
average temperature of 1659° F (1177 K). Since the other three designs had approxi- 
mately the same centrifugal stresses (29 000 to 31 000 psi or 19 995 to 21 374 N/cm^) 
and average hub temperatures (approximately 1570° F (1127 K)), their creep lives 
should depend largely on the maximum chordwise temperature differences. 

As can be seen by comparing figures 6, 7, and 9, the strut insert design has about 
a 10 percent higher temperature difference and the film- convection design has about a 
60 percent higher temperature difference than the multiple small hole design at the hub. 
Therefore, one would ejqject the multiple small hole design to be the best of the three 
and the film convection design the worst as is indeed the case. 

As expected, the creep lives based on the initial conditions were too pessimistic 
while those based on the average conditions were too optimistic. Neither gave the same 
ranking of designs as the stress relaxation analysis. The stress relaxation effects can- 
not be determined from a mere inspection of the Initial conditions because creep rates 
varied throughout the blade section and changed as the stresses become redistributed. 
Furthermore, the critical location where the worst creep occurs sometimes changed be- 
cause of the stress relaxation. 

The relative merit of the designs cannot be determined by neglecting the thermal 
stresses and considering only the mechanical stresses and the average temperature 
levels. Thus, on the basis of the average conditions, the film- convection cooled design 
appeared better than the multiple small hole design in table I whereas the stress relaxa- 
tion analysis showed the reverse to be true. 


Stress and Strain Distributions 

In figures 12 to 16 the initial stress distributions and the final stress distributions 
at the hub sections at the time when 1 percent creep strain or 100 000 hours was reached 
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are shown in part (a) and the final creep strain distributions are shown in part (b). 

Since the tolerance of the blade material to any stress level is a function of temperature 
the critical location is determined from the creep strain distribution in part (b) of fig- 
ures 12 to 16. 

Stru t insert design. - The stresses in the strut insert were higher than those in the 
shell and the stresses in the pressure side of the shell were higher than in the suction 
side, as shown in figure 12(a). However, significant stress relaxation occurred in these 
regions with compensating increases in the stresses on the shell leading edge and suc- 
tion side. 

Considerably more creep stram took place toward the trailing edge than the leadir^ 
edge. The 1-percent creep strain was reached at the extreme trailing edge (fig. 12(b)). 
This design as well as others that showed unbalanced creep strain distributions at the 
hub section would be improved by a more uniform wall temperature distribution. In this 
case, the creep life could be somewhat increased by letting the leading edge run hotter 
and reducing the metal temperature at the trailing edge. Sometimes this may be diffi- 
cult to do without adversely affecting the blade life at some other location. For example, 
with most of the blades considered the leading edge was critical at the midspan location 
and any further increase in temperature at that position would probably have a deleteri- 
ous effect on the low cycle fatigue life. 

Film and convectio n cooled design . - The webs were the coolest parts of the film- 
convection blade and, therefore, had the highest stresses. Some relaxation of these 
stresses is shown in figure 13(a). hi addition, the final stress distribution showed a re- 
duction of the trailing edge stresses. Even with this stress relaxation the trailing edge 
remained critical. The final creep strain distribution of figure 13(b) was radically un- 
balanced with most of the creep occurring at the trailing edge. 

In order to investigate the effect of the web temperature on the creep life, another 
version of this blade was analyzed in which the web temperatures were arbitrarily made 
the same as the temperatures at the adjoining parts of the shell. This lowered the web 
stresses but did not significantly affect those at the trailing edge which remained critical. 
The raising of the web temperatures also increased the average temperature level of the 
hub section slightly. This resulted in a slower stress relaxation at the trailing edge 
which reached 1 percent creep in 1300 hours compared to the 1750 hours for the original 
design (without restoring moments being considered). 

The creep life from the stress relaxation analysis was about 15 times that based on 
initial conditions compared to a factor of about 3 for the strut insert design in table I. 
This was due to the fact that the initial critical location for the film- convection design 
was at a highly stressed web node which relaxed rapidly and moved the final critical lo- 
cation to the trailing edge, hi contrast, the critical region of the strut insert design re- 
mained around the trailing edge which showed little relaxation in figure 12(a). 
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ey Stress Creep strain, 

psi (N/cm^) in. /in. 

0 Under 15 000 (Under 10 342) Under 0.001 

1 15 000 to 20 000 (10 342 to 13 790) 0.001 to 0.002 

2 20 000to 25 000 (13 790 to 17 237) .002to0.003 

3 25 000 to 30 000 (17 237 to 20 684) .003 to 0.004 

4 30 000to 35 000 (20 684 to 24 132) .004to0.005 


Key Stress Creep strain, 

psi (N/cm^) in. /in. 

5 35 000 to 40 000 (24 132 to 27 579) 0.005 to 0.006 

6 40 000 to 45 000 (27 579 to 31 026) .006 to 0.007 

7 45 000 to 50 000 (31 026 to 34 474) .007 to 0.008 

8 50 000 to 55 000 (34 474 to 37 921) .008 to 0.009 

9 55 000 to 60 000 (37 921 to 41 368) .009 to 0.010 

s Above 0.010 



(b) Creep strain distribution. 

Figure 16. - Stress analysis results for liquid-metal- and convection-cooled design. 


Transpiration cooled design . - The trailing edge of the strut developed the highest 
creep strain. This was despite the sharp stress relaxation at the trailing edge projec- 
tion shown in figure 14(a). Aside from this region, the creep strains in the strut ap- 
peared fairly well balanced (fig. 14(b)). Transpiration cooling requires a porous mesh 
material which will be resistant to oxidation at the 1600° F (1144 K) metal temperature. 
Otherwise the stress relaxation results for this design would be rather academic since 
oxidation would progressively close the fine pores and change the metal temperatures 
during steady-state operation. 

Multip le smal l hole des^n. - This configuration had the highest creep life next to 
the transpiration cooled design. Of all the des^ns it had the best strain balance be- 
tween leading and trailing edges as shown in figure 15(b) and is the closest to being 
optimized from a creep rupture standpoint. This balanced strain distribution was a 
result of the effect of creep in relaxing the leading edge and web stresses and raising 
the trailing edge stresses (fig. 15(a)). These relaxation effects are so complicated that 
the only way to deliberately design the blade with an optimum coolant flow distribution 
to give the maximum creep life would be by a trial and error procedure using a stress 
relaxation analysis. However, under some conditions where the creep life is consider- 
ably more than required for the design life of the blade, it may be desirable to redistri- 
bute the coolant to lower the temperatures of other potential problem areas such as the 
leading edge in the midspan region which may have a low cycle fatigue problem. 

Liquid met a l and convection cooled design . - The highest stresses generally 
occurred adjacent to the coolant feed channels near the leading and trailing edges. Fig- 
ure 16(a) shows only a relatively small relaxation of these stresses. The creep strain 
distribution in figure 16(b) is badly balanced with the trailing edge showing virtually no 
creep. This design could be improved somewhat by using less air to cool the trailing 
edge and more to improve the heat exchanger efficiency. 


Multiple Small Hole Blade Design Curves 

Curves of percent coolant to total flow against average blade metal temperatures 
for the multiple small hole configuration are shown in figure 17 for turbine inlet gas 
temperatures of 2270° F (1517 K) and 2570° F (1683 K) with a blade coolant inlet tem- 
perature of 1270° F (961 K) which was used in references 1 to 3 for the previous des^n 
cases. Curves are also shown for a 2570° F (1683 K) turbine inlet gas temperature with 
1000° F (811 K) and 600° F (589 K) coolant blade inlet temperatures. These curves were 
established from analytical temperature data given in reference 2 and, where these 
were not available, by using the transpiration cooling effectiveness parameter derived 
in reference 2 even though this is a film cooling des^n. It was found that this param- 
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Figure 17. - Effect of turbine gas inlet and coolant temperatures on metal temper- 
ature for multiple small hole design. Cruise condition: altitude, 75 000 feet 
( km); speed, Mach 3. 

eter gave a better correlatioji with the analytical data shown in figure 17 than the more 
frequently used convective cooling effectiveness parameter. Metal temperature for 
other conditions can be determined in a similar manner. 

Curves of creep lives as a function of average metal temperatures are presented 
in f^re 18 for centrifugal stress levels of 20 000, 24 000, and 30 000 psi (13 790, 

p 

16 547, and 20 684 N/cm^) for the multiple small hole design. Results for other centri- 
fugal stresses can be obtained by cross plotting the curves. 

These curves were determined by (1) raising all the nodal temperatures of fig- 
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Figure 18. - Effect of metal temperature on creep life of multiple small hole design 



ure 9 (p. 18) by an amount equal to the difference between the new blade average metal 
temperature and the previous one of 1629° F, (2) applying the desired centrifugal stress 
level and recomputing the creep life from the stress relaxation analysis for one refer- 
ence case, (3) calculating an effective Larson- Miller parameter, Tj^(25. 56 + logj^Qt), 
based on this case, and (4) using the effective parameter to calculate the time to failure 
at any other metal temperature, hi this way, the theoretical curves of figure 18 were 
established. The creep lives for a number of other cases at various centrifugal 
stresses and metal temperature levels were computed directly from the stress relaxa- 
tion analysis and are shown by data points in figure 18. As can be seen agreement with 
the theoretical curve was excellent, hi all the cases recomputed for other temperature 
levels, the final creep strain distribution remained essentially unchanged. The design 
curves of figure 18 do not include the effect of the centrifugal force restoring moments 
due to bending. 

As an example of how the design curves can be used, the conditions required to 
attain a 1000 hour life at a turbine gas inlet temperature of 2570° F (1683 K) with a 
30 000 psi (20 684 N/cm ) centrifugal stress will be examined. From figure 18 it is 
foimd that an average metal temperature of about 1700° F (1200 K) is required for 
1000 hours life at 30 000 psi (20 684 N/cm^). hi figure 17 it is seen that combinations 
of coolant to total flow and coolant temperatures which will give 1700° F (1200 K) are 
3. 45 percent of 1270° F (961 K) air, 2.18 percent of 1000° F (811 K) air, and 1. 42 per- 
cent of 600° F (589 K) air. Another example is to consider the effect on blade life of 
increasing the turbine inlet temperature from 2270° to 2570° F (1517 to 1683 K) while 
using 2. 8 percent coolant at 1270° F (961 K) with a 30 000 psi (20 684 N/cm^) centrifugal 
stress. This increases the average metal temperature from 1615° to 1761° F 
(1153 to 1234 K) (fig. 17) and reduces the life from over 10 000 to 140 hours (fig. 18). 

If the centrifugal stress level for the 2570° F (1683 K) condition is 20 000 psi 

(13 790 N/cm^) instead of 30 000 psi(20 684 N/cm^), the life increases to 3800 hours or 

about 27 times for a 33 percent reduction in stress level. 

As a check on these design curves a 2570° F (1683 K) turbine inlet temperature 
case with a 30 000 psi (20 684 N/cm ) centrifugal stress was analyzed using the temper- 
ature distribution given in reference 2 for a 2. 8 percent coolant to total gas flow. The 
analytical creep life was 70 hours compared to 138 hours that would be predicted on the 
basis of figures 17 and 18. The shorter analytical life was due mainly to a final imbal- 
anced creep strain distribution which had most of the creep occurring at the trailing 
edge and little at the leading edge. Another factor which serves to make the results in 
figure 18 somewhat unconservative is that the analytical temperature distributions of 
reference 2 show a tendency for the temperature difference to increase with increasing 
turbine inlet temperature. However, little effect on the temperature difference is shown 
in reference 2 when the turbine inlet temperature is kept constant but the coolant flow 
rate is chained. 
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CONCLUDING REMARKS 

This analysis does not take thermal fatigue into account and, therefore, does not 
consider the effects of discontinuities such as film cooling holes and fin to shell junc- 
tions in reducii^ the cyclic lives of the blades. The Influence of these discontinuities 
depends on factors such as their geometries, the thermal and mechanical stress fields, 
coupling due to strain gradients, and stress relieving effects due to arrays of openings 
or intersecting passages and to plastic flow. Consideration of these effects is beyond 
the present state of the art of turbine blade stress analysis except by en^irical methods 
based on extensive experimental results. The transpiration cooled blade which showed 
the longest creep life requires a porous material which could provide the fine passages 
required for ideal transpiration cooling and still tolerate thousands of hours of opera- 
tion at 1600° F (1144 K). Otherwise, long before the theoretical creep life could be 
reached the porous skin would undergo progressive oxidation which would clog passages, 
restrict coolant flow, increase metal temperatures, and accelerate further oxidation. 

Up to the present time transpiration cooling 1000 hours of steady state operation has not 
been feasible for metal temperatures much beyond 1600° F (1144 K). However, 
development of more oxidation resistant materials or coatings for porous shells could 
enable this cooling method to be used for much higher skin temperatures. Another 
problem that could limit the lives of strut supported transpiration cooled blades is the 
difficulty of obtaining a strong reliable joint at the shell to strut interface. Weak joints 
would be susceptible to thermal stress fatigue. 

The multiple small hole blade showed the best creep life next to the transpiration 
cooled blade and would appear to be the most promising of the designs for lor^ time 
steady-state operation with high metal temperature levels at the present time. The 
main reason for its superiority over the strut insert and film-convection blades was the 
more uniform temperature distribution at the hub section. However, the results of ref- 
erences 1 to 3 demonstrated that all three of these designs had fairly close overall tem- 
perature levels and temperature differences over the airfoil surfaces. Thus, for the 
2. 6 percent coolant flow ratio the range -of maximum temperature differences over the 
shell was 245° to 301° F (136 to 167 K) and the average shell temperatures ranged from 
1629° to 1653° F (1160 to 1174 K) for these designs. The steady-state lives depend on 
the mechanical loads, temperature levels, and temperature gradients. It therefore 
appears possible that, neglecting fabrication problems, these tlu*ee blade configurations 
could be designed to possess approximately the same steady-state operating lives if the 
wall temperature distributions were optimized. Therefore, the creep life for the 
multiple small-hole blade may indicate what the strut insert and film-convection cooled 
blade also attain. However, the 4 or 5 to 1 ratio between the creep lives of the multiple 
small hole and film- convection blades, despite a lower centrifugal stress assigned to 
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the latter, demonstrates the drastic effect on blade life that can occur as a result of 
differences in metal temperature distributions. 

Since the results of the heat transfer analyses of references 1 to 3 show, only a small 
effect of configuration geometry for the three basic configurations cited above it would 
seem that, as loi^ as the coolant can be effectively dispersed around the blade, the selec- 
tion of a configuration should depend largely on considerations of structure, fabrication, 
aerodynamics, and, perhaps, cyclic life. From structural and fabrication standpoints 
the strut insert design may not be desirable because of the horizontal fins which are dead 
we^ht and difficult to cast. Since both the multiple small hole and film and convection 
cooled blades possess arrays of holes near both leading and trailing edges there may be 
little to choose between them from aerodynamic and cyclic life standpoints. 

The liquid metal design had the lowest creep life because of the high metal temper- 
ature level at the hub section. This configuration had a higher overall blade temperature 
of approximately 100° F (56 K) compared to the multiple small hole, strut insert, and 
film-convection cooled configurations for the same amount of compressor discharge air. 
Therefore, this conf^ration could not be improved qualitatively by improving the tem- 
perature distribution. The higher metal temperatures probably result from added losses 
by not using the compressor discharge air to cool the blade directly. Instead, most of 
the air was used to cool the heat exchanger tubes and condense the potassium. A better 
heat exchanger design would improve the blade life. However, this design has other ob- 
vious problems. The centrifugal loads of the relatively heavy blades together with the 
large amount of material which has to be removed at the rotor rim for the heat exchang- 
ers would require a thicker and, therefore, heavier rotor. At first inspection the liquid 
metal concept does not appear too appealing. 

A reasonable goal for the operating lives of turbine blades for advanced airbreath- 
ing engines would be 1000 hours. All of the blade designs except the liquid metal config- 
uration met this requirement for the Mach 3, 75 000 feet (22 859 km), 2270° F (1517 K) 
turbine gas inlet temperature engine cruise conditions. The multiple small hole blade, 
in particular, had such excessive creep life for these conditions that it might be benefi- 
cial to decrease it by letting the leading edge run colder and the trailing edge hotter. 

This would probably help the fatigue life of the leading edge region which could be a 
more life limiting factor than the steady-state operation. 

From the design curves of figures 17 and 18 about 3. 5 percent coolant to total gas 
flow would be required for the multiple small hole blade to have 1000 hours life at a 
2570° F (1683 K) turbine gas inlet temperature under cruise conditions with a 30 000 psi 
(20 684 N/cm ) centrifugal stress. It might seem that this blade life could be increased 
by doing more cooling near the blade base and less above the midspan region. However, 
the anal 3 dical results of references 1 to 3 show metal temperatures above 1800° F 
(1255 K) and sometimes as high as 1900° F (1311 K) in the midspan to tip regions. 
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These temperatures are already dangerous from both structural and oxidation consider- 
ations and cannot be raised further without compromising the blade integrity. 

K should be emphasized that the design curves of figures 18 are an indication only of 
the blade potential. They presuppose a coolant distribution which will give an evenly 
balanced creep strain distribution between leading and trailing edge. From an analytical 
standpoint this would be difficult to achieve since the stress relaxation process is ex- 
tremely complex and it is not obvious what initial stress distribution to design for. Even 
if this could be determined by a trial and error procedure using a stress relaxation pro- 
gram, problems involved in fabrication, oxidation effects during operation, variations in 
gas temperature profile with combustor location, errors in determining blade gas pres- 
sure distributions, etc., could make this impossible to attain in practice. It, therefore, 
is probably necessary to either design for a safety factor on life or to design for a de- 
liberately unbalanced creep strain distribution. 


CONCLUSIONS 

The following conclusions were drawn from the analytical investigation of the steady- 
state operating creep lives of five turbine blade cooling designs for advanced airbreath- 
ing engines: 

1. The longest creep life was shown by the transpiration cooled design neglecting 
oxidation problems with the porous shell and a possible thermal fatigue problem at the 
shell to strut joint. 

2. The multiple small hole design showed significantly longer life than the strut in- 
sert and film-convection cooled designs. The temperature distribution across the hub 
section was the main factor in determining the ranking of these three blade designs. All 
these blades were relatively close in centrifugal stresses, average temperature levels, 
and overall blade maximum temperature differences. It therefore appears that the dif- 
ferences in temperature distribution have drastic effects and that it may be possible to 
design the strut insert and film and convection cooled blades to have about the same 
creep lives as the multiple small hole blade if their temperature distributions were op- 
timized for the engine design conditions. 

3. The liquid metal blade des^n had the worst creep life of the designs which were 
studied. The bulk and complexity of the heat exchanger did not result in any advantage 
in life. 

4. All the blade configurations but the liquid metal design showed steady-state lives 
in excess of 1000 hours for the Mach 3, 75 000 feet (22. 859 km), 2270° F (1517 K) tur- 
bine inlet temperature advanced airbreathing engine design conditions which were con- 
sidered. The multiple small hole design also had the capability of operating at 2570° F 
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(1683 K) turbine inlet temperature with a small rise in coolant flow. 

5. Raismg the web temperatures of the film and convection cooled blade to those of 
the adjoining parts of the shell only lowered the creep life. This was because the aver- 
age temperature level at the hub section increased and caused a slower stress relaxation 
to occur at the critical trailing edge. 

6. The creep lives of the configurations cannot be determined either qualitatively or 
quantitatively from the initial or the average conditions of stress and temperature. 

Stress relaxation causes the higher thermal stresses to decrease, changes the stress 
distribution and, therefore, the creep rates, and sometimes causes a shift in the critical 
region of the airfoil section. 

7. Failure to consider the restoring moments due to centrifugal force as the blades 
bend because of creep can result in calculated lives which are overly conservative. The 
application of this correction to the stress relaxation analysis caused an increase in life 
for all the blade designs although by varying amounts. The transpiration cooled blade 
life increased from 47 900 hours to an infinite life. On the other hand, the liquid metal 
blade showed only a minor increase in life. The qualitative ranking of the designs was 
not affected by the correction. 

8. The theoretical effect of an increase in the average metal temperature level can 
be predicted using the Larson- Miller parameter if the life at only one average tempera- 
ture level is known, provided neither the centrifugal stress nor the temperature distri- 
bution is changed. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 10, 1969, 

720-03-00-74-22. 
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APPENDIX A 


STRESS RELAXATION ANALYSIS 

The usual assumption of strain compatability is made that plane sections remain 
plane. For a blade with three degrees of freedom, the total normal strain is defined by 

= Cj + Cgy + CgX (Al) 

where Cj, C 2 , and Cg are the intercepts and slopes of the plane. 

The total strain on any element consists of (1) the elastic strain 
thermal deformation (3) the residual strain from previous plastic flow or the 

plastic component of strain 6p|j, and (4) the total accumulated creep strain e^j^. There- 
fore, the total strain at any time is 


^tk ■ 



“k'^k + ^pk + ^ck 


(A2) 


or 


^zk ®k^^tk “ “k'^k ~ ^pk “ ^ck^ 

Substitution of (Al) into (A3) gives the compatibility equation 

‘^zk " ®k^^l + ^2^k ^3^ ■ "k'^k " ^pk " ^ck^ 

The boundary conditions are 

= F (A5) 

yyazkXkdxdy = My (A7) 


The applied moments and My consist of components with respect to the co- 
ordinate axes of moments due to both gas and centrift^al forces. These moments are 
positive when they cause an increasing stress with increasing values of x or y. 
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Substitution of (A4) into (A5) to (A7) gives equations (A8) to (AlO), respectively 


(Cj + Cgy + CgX - ofT - £p - e^)E dx dy = F (A8) 

ff (Cj + Cgy + CgX - aT - - e^)Ey dx dy = (A9) 

ff (C^ + Cgy + CgX - ojT - €p - €^)Ex dx dy = My (AlO) 

or equations (All) to (A13), respectively. 

Cj yy* E dA + Cg yy Ey dA + Cg yy ex dA = F + yy (q-T + ep + €^)E OA 

(All) 


yy Ey dA + Cg yy Ey^dA + Cg yy Exy dA = + yy (oT + £p + e^.)Ey dA 

(A12) 

c yy Ex dA + Cg yy ^xy dA + cg yy Ex^dA = My + yy (ofT + ^ 

(A13) 


These equations can be further simplified if the coordinate axes are converted to 

centroidal axes where the centroid is defined by y Ex dA = y Ey dA = 0 and the axes 

are defined by y Exy dA = 0. The equations for convertii^ from the arbitrary set of 
axes to centroidal axes are given by equations (A14) to (A18) where equations (A14) to 
(A16) are based on reference 6. 


X 


yy ExdA 

y/EdA 


y 


yyEy dA 

yyEdA 


(A14) 

(A15) 
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(A16) 


e = -tan"^ 


^(yj ^ dA - xy yy E dA^ 


[ yy Ex^dA - ^yy e dA - yy Ey^dA + yy e dA 


xj^ = (Xjj. - x)cos 0 + (yjj - y)sin 0 


(A17) 


y{j = (y^ - y)cos 0 - (3^ - X)sin 0 


(A18) 


K Xj^ and yj^ are replaced by x|^ and y{^ the coefficients C^, C 2 , and Cg can 
be solved from equations (A19) to (A21). 


Ci = 


F + y y (ofT + €p + e^)E dA 


y/EdA 


(A19) 


^2 = 


X +yyioi'^ + ^ 
yy E(y’)^dA 


(A20) 


_ ^ + yy ^p + ^ 

yy E(x*)^dA 


(A21) 


Substitution of equations (A19) to (A21) into (A4) gives which in summation 
form is 
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n 


^ Z (“A * «pk * Sk)®k 

k=l 


n 


Z =k^i 


k 


k=l 


+ y’Eb 


n 


^ ^“k'^k ^pk ^ck^^k^k ^k 


k=l 


n 


Z ^ 

k=l 


+ x’E. 


n 


M 


y E (“k'^k + ^pk + ^ck)®k^k ^1 


k=l 


n 


S ^k 

k=l 


- ®k<“kTk + 'pk + "ck> <'^22) 


The strain hardening rule requires that the change in creep strain over a time in- 
terval At, be determined at a constant stress which is calculated from the previous ac- 
cumulated creep strain. Stresses are calculated from equation (A22) and creep strains 
from equation (1), The procedure is as follows: 

(1) Calculate from equation (A22). This process is started at t = 0 hours 
where = 0. 

(2) Confute a fictitious value of time t’ for each node by solving equation (1) for 
t’ using from step 1 with the same used in step 1. 

(3) Solve for the new creep strain at the end of the next time interval by sub- 
stituting t’ + At for t in equation (1), 

(4) Go back to step 1 using the new 

This procedure is repeated over successive time intervals until a creep strain limit 
is reached. The total elapsed time is the sum of all the time increments which were 
applied. Small time increments are used in order to apply an incremental method. 
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APPENDIX B 


ANALYSIS OF RESTORING MOMENTS DUE TO CREEP 

The blade centrifugal force distribution in the spanwise direction was approximated 


by 


F = F 
z z 




(Bl) 


The curvature at any span location is assumed to be defined by 

R R \L / 
o 

The usual approximation for curvature is made that 


1 _ dV 


The second and first derivatives of the deflection along the blade span are then 
given by 


(B2) 


(B3) 


d^w 


dz2 

R 

( 

dw _ 

L 

dz 

1^0- 






(B4) 


(B5) 


The restoring moment due to the centrifugal force is 

M = - J''P^ dw 

Substitution of equations (Bl) and (B5) into (B6) yields 


(B6) 
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200(|-jj^ + 0. - 1. L 


which results in 


M = -0. 118 


F L“ 
z.o 


Since l/R^ is equivalent to the slope of the hub section, the restoring moments 


become 




My = -0.118F^^^L"C3 
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